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Maximum Likelihood Estimates of Position Derived From 
Measurements Performed by Hyperbolic Instruments^ 

Eugene Lukacs 

Certain electronic surveying instruments — referred to as hyperbolic instruments — 
determine the difference between the distances from two fixed stations to a moving object. 
Data obtained from a single observation yield one relation between the three coordinates of 
the moving object and restrict it to a hyperboloid of revolution. To locate a position in 
space, observations from at least three pairs of stations are needed. If more than three 
observations are taken, more than three equations for the three coordinates are obtained. 
In general these equations will not be compatible; it is then a statistical problem to find an 
estimate for the unknown position. In this paper maximum likelihood estimates are 
obtained for positions in space derived from such hyperbolic data. 



1. Introduction 

Several ehn'tronic surveying instruments have 
been developed in i-eeent years. These instruments 
utilize the Doppler effect ^ ^ and are (h^sio-ned to 
locate the position of a moving oV)j(H*t. In this 
paper we are considering measurements ohtained 
by using certain types of these instruments, which 
are often referred to as hyperbolic instruments. 
Hyperbolic instruments, such as Raydist (see foot- 
note 2), determine the difference between the 
distances from two stations to a moving object. 
Data obtainable from one observation (using on(» 
pair of stations) yield one relation among the three 
coordinates of the moving objcM't. To locate a 
position in space from hyperbolic data, at least 
three observations are needed. If more than three 
observations are taken, more than three equations 
for the three coordinates are obtained. In general, 
these equations will not be compatible; it is then a 
statistical problem to find an estimate for the 
unknown position. The use of exactly three obser- 
vations is justified if it can be assumed that the 
data are not affected by any observational errors. 
In general, this assumption is not realistic and should 
*not be made without a thorough investigation. It 
appears desirable, therefore, to use more than the 
minimum number of observations and to estimate 
the location of the moving object by statistical 
methods. The purpose of this paper is to obtain 
maximum likelihood estimates for positions in space 
derived from hyperbolic data. 

2. Formulation of the Assumptions and the 
Problem 

We denote by A^>4 the num.ber of stations 
Sa{a=l,2, . . . N) and by (aai,«a2,«a3) the known 
coordinates of the station *§«. Let P be the unknown 
position of the m.oving object and (^i, 62, 6^) be the 

1 An unclassified treatment of the mathematical portion of a classified Navord- 
report written by the author while a member of the staff of the U. S. Naval 
Ordnance Test Station, Inyokern, Calif. 

2 J. F. McAllister, Measuring velocity of V-2 rockets by Doppler effect, 
Tele-Tech. 6, No. 2, p. 56 to 59 and 129 (1947). 

3 Charles E. Hastings, Raydist— a radio navigation and tracking system, 
Tele-Tech. 6, No. 6, p. 30 to 33 and 100 to 103 (June 1947). 



coordinates of P. A single observation based on 
the stations S« and S^ gives then a determination of 
the distance PSa—PS^. If A^ stations are available 
it is possible to make at most 



n- 



_N(N-1) 



observations based on different pairs of stations. 
However, not all these observations will be used in 
practice. Let r^^ [A^(A^— 1)]/2 be the nimiber of 
observations actually taken. We make the follow- 
ing assumptions: 

(I) The hyperbolic m.easurements Za^ are observa- 
tions on n independent random variables Za0. The 
variates Z^^ have rectangular distributions with mean 
11 a0 and comm.on range 2r. Further Z^a^—Za^ and 

(II) The parameters of the n rectangular distribu- 
tions m.ay be expressed in terms of the coordinates 
of the unknown point by m.eans of the formula 



(III) The distances da and d^ from the position of 
the moving object at time t=0 to the stations Sa 
and S^, respectively, satisfy the inequality 



da — d^^r. 



(2) 



The nature of the problem., as well as its solution 
are to a great extent determined by assumption (I). 
To justify this assumption it is necessary to describe 
briefly how hyperbolic data are obtained. 

It is desired to determine the instantaneous loca- 
tion of a moving object. A transmitter is carried 
in the object and emits a continuous signal of fre- 
quency /o. A second transmitter (the reference 
transmitter) emits a signal of a frequency /'o? slightly 
different from/o- Each station receives both signals 
simultaneously. A certain beat frequency ria is 
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thus associated with each station Sa', Ua depends on 
the velocity component of the object away from Sa 
and on the frequencies /o and /'o- To obtain an 
observation two stations Sa and S^ have to be hnked 
by wire or by radio The frequencies Ua and n^ 
give then a beat frequency; the quantity recorded is 
the num.ber of cycles Ua^ of the beat frequency signal 
occurring during the tim.e of observation. It can be 
assum^ed that the observation Za^ on Za^ can be satis- 
factorily approxim^ated by a linear function of ria^ 
given by 



Jo 



(3) 



In this formula c is the speed of light. The frequency 
Ua^ m.ay be recorded continuously to permit reading 
Ua^ to fractions of a cycle. Alternatively it is possible 
to devise a counter that registers full cycles only.* 
In either case it is easily possible to determine a 
range for the reading error of ria^ (and therefore also 
of Za^). Moreover, the recording devices are of such 
a nature that any error within the possible range is 
equally likely to occur. It seem^s, therefore, quite 
realistic to assume that the errors have a rectangular 
distribution about zero with a certain range. This, 
however, is equivalent to assumption (I). Equation 
(2) is a condition im.posed on the location of the 
stations. As the half range r is rather sm.all (2) 
implies hardly m.ore than da^d^. This considera- 
tion shows that (2) is not a very serious restriction. 
The aim of this paper is the determination of a 
maximum likelihood estimate of the unknown posi- 
tion. That is, the unknown coordinates ^i, 02, d^ 
have to be determined so that the greatest possible 
probability is given to the actual observations. This 
discussion will be based on assumptions (I), (II), 
and (III). 



3, The Maximum Likelihood Estimate 



Let/a^(s) be the frequency function of the random 
variable Z^^. According to assumption (I) 



1 



if fJia^—r^ z ^ Ma^+r. 



,0 otherwise 

The likelihood function ^ assumes, therefore, also 
only two values, the value (l/2r)'' and the value zero. 
If one wants to determine the unknown parameters 
^1, ^2, ^3 so that the greatest possible probability is 
given to the actual observations, one has only to 
chose the parameters so that the likelihood function 
is not zero. This is the case if and only if each 
observation Za^ falls within the range of Z^^. 
That is, the likelihood condition 

* Hastings Bulletin R-23, Complete Raydist tracking system for missiles (Sept. 
1947). 

5 A description of the maximum likelihood method and a description of the 
likelihood function may be found in, H. Cramer, Mathematical Methods of 
Statistics, p. 498 (Princeton Univ. Press, Princeton, 1946). 



Ma^ — ^^^a/3^Ma^ + ^, 



(4) 



must be satisfied for the n pairs {a,^) used. From 
(1) it is seen that the likelihood condition (4) consists 
of a set of inequalities involving di, 62, 6^. A point 
(^1, 02, 6z) of the parameter space is said to be a 
maximum likelihood estimate if its coordinates 
satisfy the maximum likelihood condition. In gen- 
eral, inequalities (4) will not have a unique solution. 
Either all the points of a certain region in the para- 
meter space will satisfy (4), or it may happen that 
no maximum likelihood estimate exists. If maxi- 
mum likelihood estimates exist, one may consider 
the set of all the maximum likelihood estimates. 
This set will be called the maximum likelihood region. 
Any point in the interior of the maximum likelihood 
region may be used as an estimate for the unknown 
position. The most important problem is the deter- 
mination of the maximum likelihood region. This 
problem will be solved by expressing the boundaries 
of this region in terms of the actual observations. 
This requires certain transformations of the likeli- 
hood condition (4). If we substitute the value of 
jua^ from (1) and transpose we obtain 



^a^—r^^^{ei — aaiy—-J^{0i — a^iy^Za^+r 

(4a) 

We next introduce the following notation 



^0=—— (da — d^—r) 

n, = ^^(da-d^-\-r) 
c 



(4b) 



From (4a) we see ni = nQ—2rJo/c<Cno. 

In the following we distinguish three cases, namely, 



Case (A) 
Case (B) 
Case (C) 



na0<ni 



From (3) it is seen that 

Za^—T^O and a fortiori ^a^+r> in case (A) 
Za^+r<CO and a fortiori s^^— r<0 in case (B) 

^a^+r"^ ^ Za^—r in case (C). 



y (5) 



We first consider the conditions imposed by a single 
observation Za^ based on an arbitrary pair of stations. 
It is no loss in generality if we assume a=l and /S=2. 
To simplify the discussion we choose an appropriate 
system of coordinates. Let the position of the mov- 
ing object at time ^=0 be the origin and choose the 
^i-axis so as to be parallel to the line joining Si and 
S2 the ^2-axis should be in the plane OS1S2 and per- 



198 



pendicular to 8182.^ The ^3-axis is perpendicular to 
the ^1^2-plane; finally some orientation is given to 
the axes. With this choice of coordinates, we have 



(Il3 = (l23'^0 



(6) 



We shall write an = ai, an^a^; ai2=a22=b. 

To make the formulae more concise, we introduce 
the following notation: 



2 = |si2-r|; 
2 = \zi2+r\; 



r=\Z2 



-r\, 



2 =|221+?-| 



(7) 
(7a) 



i?i=PS(=V(ei-a,r+(e2- 6)2+0^ for i = l,2, (7b) 



wi = (ai+«2)/2, 
d = {ai—a^j2; d* = (a2—ai)l2, 



T.= m-^; 



T,=m-i-T 



Id' 



T* = m- 



Tt = m- 



'4d*' 
' 4d*' 



A 2fi? 



2d* 



A*=^ 



A2 = 



2d 



At= 



2d* 



(7c) 
(7d) 

(7c) 
(7f) 
(7g) 
(7h) 



Considering that according to assumption (I) 
Z2i=—Zn, we obtain from (7) and (7a) 



£*=z and 2*=z. 


(8) 


We see also 




d*=-d, 


(8a) 


and therefore 




Tt + T2=T,+n=2m, 


(8b) 


from these relations we see 




A^= — A2 and At= — Ai. 


(8.-) 



A simple computation gives Bi^—R2^=—4:d(6i — m) 
so that we have with the aid of (7e) and (7f) 

6 If the origin is located on the line S1S2 the three points OS1S2 do not determine 
a plane. In this case the ^2-axis may be chosen in an arbitrary plane so as to 
intersect the 0i-axis perpendicularly. 



m-Rl-i'=4:d(T2-e,) 
R\-m~l'' = U{Ti-d,) 
Rl-Rl-z_^=4:d{Bi-T*) 
Rl-Rl-T'=id{ei~Tt)J 



(9) 



We next rewrite the maximum Hki^liliood condition 
for each of the three cases (A), (13) and (C). 
Case (A) 

^ From (5), (7) and (7a) it is seen that 2 = ^12— r, 
^ = 2i2+r.^ It follows then from (4a) that 0<.2 + 
^2^^i^^+^2. By squaring and transposing we 
obtain the inequalities 



R2^ 



Rl-Rl-z' 

2z 



and R2> 



R\-Rl-1'' 

27 



Considering (9), (7g), and (7h), we obtain from these 
relations the final form for the maximum ]ik(*liliood 
conditions 

(A.l) i?2^^2(r2-^i) and (A.2) R2'^A,{T,-e,). 

In an analogous manner we obtain the likelihood 
conditions in Case (B). 

(B.l) R,^At{Tt-e,) and (B.2) Ri^ At{Tt-d,). 

It is seen that the inequalities (B.l) and (B.2) may 
be obtained from (A.l) and (A.2) by writing Ri in- 
stead of R2 and by attaching an asterisk to the A 
and T symbols. 

Case ((7) Similarly, we obtain the likelihood condi- 
tions 

(C.l)R2^A(l\-d,) and (C.2)R,^Ar{Tt-d:). 

It is seen that (C.1)=(A.2) and (C.2) = (B.2) 

4. Discussion of the Likelihood Condition 

Before discussing the likelihood condition, it is 
desirable to derive a simple relation. 

From (2) we have c?i>(i2; in our coordinate system 
this means 

aiyalor |ai|>|a2|. 
From this we see easily 

sign G? = sign (ai — a2)=sign ai = — sign c?*. (10) 
Here sign x stands for "the sign of x'\ that is, 
sign x = -\-l if x>0 and sign x= — l if x<0. 

We first discuss the case (A) and start with condi- 
tion (A.l). It is seen immediately that (A.l) can 
only be satisfied if .42(^2-^1)^ 0. 

Considering (7h) and (10) this necessary condition 
can be written 



T2 sign di^ 61 sign ai. 



(11) 



This condition restricts the point P to a half space 
determined by the observations and the location of 
the stations. In the following we assume that ( 1 1) is 
satisfied. We can, therefore, square (A.l) and obtain 



957570—51- 



-6 
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it follows that 

(e2-by+ei^{At-m-2{AiT,-a,)e^+(AiTi~ai). 

(12) 

We have to distinguish two possibihties (a) If 
Al=l then from (7h) and (7f) T2=m — d = a2. 
The right-hand side of (12) is then zero so that (12) can 
only be satisfied for the points of the line 62= b, 

We consider next the case where (b) Aly^l. From 
(7f) and (7h) it is seen easily that 



AlT2-a2=m{Al-l). 



(12a) 



Using this formula we obtain by a simple computa- 
tion from (12) 

{Ai-\){e,-my-{d,-hy-0i^^^^- (12b) 

The equality sign in (12b) determines a quadric 
surface of revolution Qi. The center of this quadric 
is the point (m,6,0); the axis of revolution is parallel 
to the ^i-axis and is therefore the line joining the two 
stations Si and ^^2. The surface Qi is a hyperboloid 
if ^2>1 and an ellipsoid if Al<i\. Denote by u 
the real (resp. major) semiaxis and by v the imaginary 
(resp. minor) semiaxis of ^1. Then it is seen from 
(12b), (7h), and (7f) that 



u'^=^ and v^=^ Al— 1 . 
4 4 ' ^ ' 



(13) 



Let 2e be the distance between the two focal 
points of Qi, then 6^=='?^^4-^^sign(^2— 1) or from (13), 
e=d and u=3/2. 



e=d and u=— z. 



(13^a) 



The focal points of the quadric Qi are therefore the 
stations Si and S2. The quadric surface Qi divides 
the space into two regions, we call the region con- 
taining the center the interior of ^1 and its comple- 
ment the exterior of Qi. We can then say that 
condition (A.l) is satisfied outside Qi if Af^l and 
inside Qi if Al<il. If Al^l then (A.l) is satisfied 
only on points of the line 62=}), ^3=0 for which 
A = sign(r2-6iO. , 

We next consider condition (A.2). This ine- 
quality is satisfied in two different regions. In the 
first region the right-hand side of (A.2) is positive, 
in the second the right-hand side of (A.2) is negative 
or zero. We assume first that the right-hand side 
of (A.2) is positive. By squaring and transposing 
we obtain then 

{e2-hy+el'^{A\-l)e\-2{AlTi-a2)ei + A\Tl-al 

(14) 



Again we distinguish two possibilities: (a) If A\=\ 
then from (7e) and (7g) Ti=m-d=a2. The right- 
hand side of (14) is then zero; in this case (14) is 
satisfied for all points; (b) if AIt^I we see easily 
from (7e) and (7g) that 



AlTi-a2=m{Al-l). 



(14a) 



Using this formula (14) reduces by a simple compu- 
tation to 

{A\-\m-my-{e,-hy-gi^^^^=^- (i4b) 

The equality sign in (14b) determines again a 
quadric surface of revolution Q2. The center of this 
quadric is the point (m, 6, 0); the axis of revolution 
is again the line joining the stations Si and S2. The 
surface Q2 is a hyperboloid if ^i^> 1 and an ellipsoid 
if Ai^<i\. Denote by p the real (major) semiaxis 
and by ^ the imaginary (minor) semiaxis of (?2. 
Then it is seen from (14b), (7g) and (7e) that 



f=^ and g_^==— 



Al-\l 



(15) 



Let 2/ be the distance between the two focal points 
of ^2. By an argument similar to the one used 
before, we obtain 



f=d and ^=7 



(15a) 



so that the quadrics Qi and Q2 are confocal. We 
summarise these results in the following statement. 
If the right-hand side of (A.2) is positive, the 
condition (A.2) is satisfied inside the quadric Q2 
when Ai^y>\, but outside ^2 when ^i^<l. In 
case Ai^=l, condition (A.2) is satisfied for all points 
for which A = sign (Ti — di). 
We investigate next the second region in which 
(A.2) is valid by assuming that the right-hand side 
of (A.2) is nonpositive; that is 



Ai(Ti- 61)^0. 



(15b) 



From (7g) and (10) we see that sign Ai=sign ai so 
that (A.2) is satisfied if 



Ti sign ai^ ^i sign ai. 



(15c) 



From (15c) and the necessary condition (11) it is 
seen that (A.2) is satisfied if di belongs to the interval 
(Ti, T2) in case (15b) holds. Therefore, (A.2) is 
satisfied either in the above-mentioned region deter- 
mined by Q2 respectively in a half space or in a region 
bounded by two planes perpendicular to the line 
joining Si and S2. 

In case (A) we have ^=Zi2—ry>0 so that z>f. 
Therefore, we see from (7g) and (7h) that A2^>Ai^ 
and using also (11a) and (15a), u<^p. We are now 
ready to determine the region where (A.l) and (A.2) 
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are satisfied simultaneously. We have to distiiiguisli 
five possibiUties. 

1. A2^>Ai^>l. It is easily seen that the point 
iT2,b,0) [resp. the point (Ti,b,0] is on the line seg- 
ment joining the vertices of Qi (resp. Q2). From 
our previous results we see that (A.l) and (A.2) are 
satisfied in one-half of the hyperbolic shell bounded 
by Qi and ^2. 

2. A22>l = Ail The point {T,,bfi) coincides 
with S2, while (^2,6,0) is on the line segment joining 
the vertices of ^1. Conditions (A.l) and (A.2) are 
satisfied in one-half of the outside of the hvper- 
boid ^1. 

^ 3. A2'>l>Ail The point (^1,6,0) is then out- 
side ^1 and outside Q2, while {T2,b,0) is inside <?i. 
The region consists again of one-half of the outside 
of the hyperboloid Qi. 

4. A2^=l>Ail In this case conditions (A.l) and 
(A.2) are satisfied only by points of the line S1S2, 
which obey (11). 

5. If l>A2^>Ai^ no region exists where (A.l) and 
(A.2) are simultaneously satisfied. It is worth while 
to remark that cases (2), (3), (4), and (5) are of 
little practical interest. From (7g), (7h), and (10) 
it is seen that Ai^^l (respectively A2^^1) is equiv- 
alent to 

\2d\^\2i2—r\ (respectively \2d\^\2i2+r\). 

A situation where these relations hold is conceiv- 
able but unlikely as long as r is small. This is seen 
if we rem.ember that the difiference betwcHMi two 
sides of a triangle is always less than its third side. 

Case (B) can be discussed independentl}^ by repeat- 
ing the arguments used for case (A). 

We can then summarize the results by distinguish- 
ing the following five cases: 

1. Ai2>A2'>l. In this case (B.l) and (B.2) are 
satisfied in one half of the hyperbolic shell bounded 
by Qi and Q2. 

2 A 2^1^ A 2 1 conditions (B.l) and (B.2) are 

q' a^2\i\ a ^2*1 satisfied in one half of the out- 

3. Ai >1>A2 .J gj^j^ ^f ^^ 

4. Ai2=l>A2l Then (B.l) and (B.2) arc satis- 
fied only by points of the line S1S2 that obey (11). 

5. If l>Ai^>A2^ no region exists where both 
(B.l) and (B.2) are simultaneously satisfied. Cases 
(2), (3), (4), and (5) offer again little practical 
interest. 

We finally consider Case (C) with the conditions 
(C.1)^(A.2) R2^A,{T,-e,) and (C.2)^(B.2) 
i?i^^i*(ri*-0i). It is seen that (C.l) is identical 
with (A.2) and (C.2) with (B.2), therefore, the pre- 
vious results can be applied immediately. It is 
advisable to remark that no necessary conditions 
[similar to (11)] exist ; however, (C.l) as well as (C.2) 
are satisfied if the right-hand sides are nonpositive. 
Considering the previous results, case (C) may be 



summarized by distinguishing the following five cases. 

1. Al>Aiyi [or AiyAiyi]. Conditions (C.l) 
and (C.2) are satisfied in a simply connected region 
bounded by one nappe of Qi and by one nappe of Q2 
so that the convex side of one nappe is turned toward 
Si J while the convex side of the other nappe is 
turned toward S2. 

2. AiyAi=l [or Al>Al=ll The conditions 
(C.l) and (C.2) are satisfied in the region that is 
bounded by one nappe of one of the quadrics and 
that contains the center. 

3. AiyiyAl [or AiyiyAll The conditions are 
again satisfied in the region that is bounded by one 
nappe of one of the quadrics and that contains the 
center. 

4. Al=l>AnoT Al=l>Al]l^^'' conditions are 
5.l>Al>Ai[orl>Ai>4]\ everywhere satis- 

In case (C) the observation gives only little infor- 
mation about the location of the moving object. 
However, case (C) can occur only if 1^12 |<r. This 
will happen only for a small portion of the path of 
the moving object if r is small. Moreover, it should 
be possible to locate the stations in such a way that 
case (C) occurs only for one observation at a time so 
that this possibility should not cause any difficulty 
in locating the moving object. 

5. The Maximum Likelihood Region 

In section 2, it was assumed that n hyperbolic 
observations (based on n pairs of stations) were 
taken to determine the position of the moving object. 
In section 3 we derived the likelihood condition 
imposed by a single observation and discussed in 
section 4 the resulting information on the location of 
the moving object. It was shown that the moving 
object was confined to a certain region in space with 
known boundaries. We finally consider all the n 
observations. Each observation determines a region; 
if there exists a point set in the parameter space 
common to all these regions then a maximum likeli- 
hood region R exists and is identical with this point 
set. The boundaries of R are hyperboloids of revo- 
lution. Any point inside the region R can be con- 
sidered to be a maximum likelihood estimate.^ 

With an appropriate choice of stations it is easy to 
devise a computational procedure that leads always to 
an interior point of R. If the maximum likelihood 
region becomes small, then it is practically possible 
to use any point interior to R as an estimate for the 
unknown position of the moving object. 



7 The maximum likelihood region cannot increase in size when the number 
of observations increases. On the contrary, it is to be expected that it will 
decrease. 



Washington, August 18, 1950. 
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